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Abstract 

We analyze the anisotropic Kardar-Parisi-Zhang equation in general sub- 
strate dimensions d! with spatially correlated noise, (ry(k, a;)) = and 
(77(^)77(^,0/)) = 2D(k)5 d '(k + k')6(uj + io') where D(k) = D + D p k~ 2 P, us- 
ing the dynamic renormalization group (RG) method. When the signs of the 
nonlinear terms in parallel and perpendicular directions are opposite, a novel 
finite stable fixed point is found for d' < d' c = 2 + 2p within one-loop order. 
The roughening exponent a and the dynamic exponent z associated with the 
stable fixed point are determined as a = — ^-j^), and z = 2 — a. For 
d' > d' c , the RG transformations flow to the fixed point of the weak-coupling 

limit, so that the dynamic exponent becomes z = 2. 
PACS numbers: 61.50.Cj, 05.40.+j 
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Recently, there have been many studies in the field of nonequilibrium surface growth. A 
number of discrete models and continuous equations for the surface growth phenomena have 
been introduced and studied |I| . One of the most interesting features in the nonequilibrium 
surface growth is the nontrivial dynamic scaling behavior of the interface width, i.e., 

W(L, t) = ( JL J2(K - hfY 12 ~ L°f{t/L*), (1) 

i 

where hi is the height of site i on the substrate, h, L, and d' denote the mean height, system 
size, and the substrate dimension, respectively. The symbol (• ■ ■) stands for the statistical 
average. The scaling function f(x) approaches to a constant for x ^> 1, and f(x) ~ x 13 for 
i < 1 with z = a/ (3. The exponents a, (3 and z are called the roughness, the growth, and 
the dynamic exponent, respectively. 



The simplest nonlinear dynamic equation describing a growing surface with lateral growth 
effect was introduced and studied by Kardar, Parisi, and Zhang (KPZ) 0, which is given 
by 

d ^ = ^h + \{yhf + r l ^t). (2) 
The noise r/(x, t) is, in general, correlated spatially and temporally, (r)(k, u)) = 0, and 

(fj(k, u)fj(k', w')> = 2£>(k, uj)5 d '{k + k')5{u + u'), (3) 

where r?(k, u) is the Fourier transform of the noise 77 (x, t). When the noise is correlated only 
spatially, an interesting form of D(k,u) can be written as D(k) = Dq + D p k~ 2p ||. Here, 
the Dq term is needed to make the RG transformation closed. The white noise corresponds 
to the limiting case of p = 0. When p > 0, the noise is correlated in space with a power law. 
It is useful to remind that the KPZ equation is invariant under the Galilean transformation, 
which yields the scaling relation, a + z = 2 . 
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Recently the anisotropic KPZ (AKPZ) equation was introduced to study the surface 
growth on the vicinal substrate Besides, the AKPZ equation can also be applied to 
various physical growth problems such as the ion-sputtered surface growth ||, and the 
surface growth on the reconstructed surface structure [0, etc. The AKPZ equation is written 

as 

d -^> = v ± Vlh + in V\h + ^(V^) 2 + ^(V,,/.) 2 + tKx, t), (4) 

where Vj_(V||) is the gradient along the perpendicular (parallel) directions. The anisotropy 
means r„ = v\\/v±_ ^ 1 and r A = Am/ Aj_ ^ 1. Here, we consider only the case of positive u± 
and v\\ for stable surfaces. The AKPZ equation was studied by Wolf for the case of white 
noise in 2+1 dimensions using the dynamic renormalization group (RG) method He 
found that when the signs of A's are opposite, the nonlinear terms turn out to be irrele- 
vant under the RG transformation. As a result, the AKPZ equation with opposite signs of 
A's belongs to the weak coupling limit, the Edwards- Wilkinson (EW) universality class J7|. 
A stochastic model associated with this weak coupling behavior of the AKPZ equation in 
2+1 dimensions was introduced by the present authors as a generalization of the Toom 
model U . By the simulations, we confirmed that the height fluctuations width increases log- 
arithmically both with time before saturation, and with system size after saturation, which 
is the signature of the EW universality class. 



In this paper, we extend the study of the AKPZ equation to the case of spatially corre- 
lated noise, D{k) = D + D p k~ 2p , in general substrate dimensions d', using the dynamic RG 
method. Thus the previous study by Wolf corresponds to the limiting case, p = and d' = 2, 
of the current study. By this extension, we find that when the signs of the nonlinear coupling 
coefficients are opposite, a new "stable" fixed point of the RG flow exists within the one-loop 
order in the "off-axis region" of the parameter space (Uq, U p ) provided d! < d' c = 2 + 2p. 
Here U and U p are dimensionless effective coupling constants, U = K d i_ 1 D X']_/h'^_ and 
Up = Kd'-\D p X 2 L /i'\ 1 where = Sd' /(2n) d ' , with Sd' the surface area of a unit sphere in 
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d! dimensions. Moreover, the dynamic exponent z associated with the fixed point is deter- 
mined exactly as z = 2 — |(p — ^-^). Thus, even at d! = 2, the exponents a and z are 
exactly obtained as z = 2 — ^ and a = y for finite p. This result would be interesting, 
because it is rare to know the dynamic exponent exactly at d! — 2 in other surface growth 
problems. As the dimension d! approaches d' c for finite p, the stable fixed point moves to the 
origin, (Uq, U*) = (0, 0), and the nonlinear terms become irrelevant for d! > d' c . Conversely, 
as the power p approaches zero at d' = 2, the fixed point shifts towards the origin, and the 
weak coupling behavior derived by Wolf is recovered. On the other hand, when the signs 
of the nonlinear terms are the same, the behavior of the RG flows is similar to the one for 
the isotropic case with the correlated noise studied by Medina et al. ||, so that we will not 
consider this case further. 

The equation we study here is the AKPZ equation, Eq. (4), in general d' dimensions, 
among which one dimension is assigned to the parallel dimension on the substrate, and the 
remaining d! — 1 dimensions to the perpendicular dimensions. The calculations of the RG 
transformation were performed by combining the two methods, one of which was introduced 
by Medina et al. to study the isotropic case with the correlated noise, and the other by Wolf 
to study the anisotropic case with the white noise. The steps of the RG transformation are 
following. First, the anisotropic exponent \ is introduced to relate the two characteristic 
length scales as £n ~ where £n and £j_ are the characteristic length scales in parallel and 
perpendicular directions on the substrate, respectively. Next, the coarse-graining transfor- 
mation is performed within the one-loop order by integrating out the fluctuations of heights 
within small length scales, which correspond to the wavevectors, e~ l ir/a < \k±\ < ir/a and 
e^ x 7r/a < \k\\\ < ft/a, (a is the lattice constant) for the parameters, u±, vn, Aj_, An, and 
the effective coupling constants, Uq, and U p . Through this coarse-graining process, the A's 
remain unchanged, which is due to the invariance of the AKPZ equation under the Galilean 
transformation, 
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(5) 



where % denotes _L or ||, and e is an infinitesimal angle of tilt. As a final step, the rescalings 
are performed as x ± — > e l x±, x\\ — > e lx x\\, h — > e la h, t — * e iz t. Under the rescalings, one 



has A i -> e'( a+z - 2 )A , and 



3 2/(l-x) 



r A . From the scale invariance, the scaling relation, 



a + z = 2, and x — 1 are obtained, provided that \± ^ and r A 7^ 0. The case of Aj_ 7^ 
and r A 7^ is much more realistic than the case that one of the A's is zero. Fortunately 
for x — 1j the recursion relations are explicitly calculable, but they could not be done for 
X 7^ 1- We find that the RG recursion relations for x — 1 are as follows. 



= v± 



z-2-(U + U p )(A 



'l- r -±)A- . 
rj d 



2C 

d'-l 
2C 



U, 



2pB 



P Al 



d'-l 



+ 2r x E] +2pr u U P [( 



B 



d'-l 



+ 



(6) 
(7) 



1 

^ = U (2 - d') + V&G + 3A - J^) + U U P (3A - + 2G) + U> p G, (8) 



dUp 
dl 



(2-d' + 2 P ) + 3(U + U p )(A - J£_) - U p -^ 



(9) 



Here A,B,C,E,F, and G are defined as A = 7(0,1,^,2), B = 7(0,l,f,2), C = 
7(0,1,^,3), E = 7(2,1,^,3), F = 7(2,l,f,2), and G = 7(0,2,^,3), respectively, 



where 



Jo (1 + ?r) 7 (l + r v y 1 ) 6 



Thus, A, 7?, C, E, F, G are functions of the variables, d' , r V) and r\. 



Before we analyze the RG recursion relations in detail, it is interesting to note that 
Eqs. (6) and (9) have the similar expressions. To make parameters unchanged, if we set 
dr u /dl = and dU p /dl = 0, we can easily get the relation, z — 2 = —(2 — d' + 2p)/3. So the 
dynamical exponent can be written as 

2 d'-2 2 d! -2 

^ = 2-g(p-^-), a= 3^-^-)' ( U ) 
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provided u± ^ and U p ^ 0. Thus, if there exists a stable fixed point at finite values of u± 
and Up, we then immediately have the expressions for z and a. 



In order to find the stable fixed point of the recursion relations, it is of help to recall the 
RG flow for the case of the white noise, p = at d' = 2 j3J. For this case, there exist the 
fixed points r* = ±r\, which were obtained by setting dr v /dl = 0. For r* = r\, the fixed 
point of the effective coupling constant, U = DX 2 /i> 3 at zero, is unstable, while it is stable 
for r* = — r\. Thus for r\ < 0, the nonlinear terms of the AKPZ equation become irrelevant 
even though they are finite, and the AKPZ equation belongs to the EW universality class. 
Following the similar steps, we have analyzed the recursion relations for the case of correlated 
noise with d' = 2 first. When d' = 2, Eq. (7) becomes explicitly as 



dr„ _ ir{r v + r\) \ (JJ | TJ Vl t / ^2 (v , _ „ ^ j_ AnTT (ir , _l 0lr , /^r_r >lr ,l_ i i 
dl 



(U + U p ){\ + ^Y(r x - Tv ) + ApU p (r x + 2r Av ^ - 2rS - r 



'121 



8(1 + ^) 2 r3 

Thus, r* = — r\ is the solution of dr u /dl = even for the case of correlated noise. Solving 
dU /dl = and dU p /dl = with r* = — r A in Eqs. (8) and (9), we then find two fixed points 
other than that of Uq — U* = which is unstable for p > anyway. The stability analysis 
shows that, of the two new fixed points, one is unstable and the other is stable for all p > 
and r\ < 0. The stable fixed point is located at 

u; = ^3 p-9 + 3^T2cT9j (13) 

v , = 4PVK / C 2 + 4C + 9-(C + 3)VC z + 2C + 9 \ (14) 



in the one-loop approximation where 



^ ~ q + /p^2" 



Thus, we indeed have a stable fixed point in the 'off-axis' region and hence the dynamic 
exponents are given by Eq. (11) with d 1 = 2. As p decreases to zero, the stable fixed point 
moves to (C/q, U*)—(0,0), and the case of white noise is recovered. As mentioned before, the 
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recursion relations for x 7^ 1 could not be derived explicitly. Thus we consider the case of 
X 7^ 1 by extrapolating the analysis for \ = 1 with the recursion relations, Eqs. (6-9) in the 
following way. Since x ^ 1 means r\ = 0, we examine the behavior of the fixed point by 
taking the limit r\ — > 0~ for fixed p. As can be seen in Eq. (13) and (14), the stable fixed 
point moves to the origin in the parameter space as r A approaches to 0~, which suggests 
that the case of x 7^ 1 belongs to the weak coupling limit. 

For general dimensions d', the explicit expressions of the recursion relations are more 
complicated. Thus, it is difficult to derive the analytic formula for r* from dr u /dl = 0. 
Instead, the stable fixed point, (r*, Uq, U*), is found numerically by iterating the RG trans- 
formation. We find that the stable fixed point exists in the off-axis region of the parameter 
space for d' < d' c = 2 + 2p, and it locates at Uq = U* = for d' > d' c . This is in accord with 
Eq. (9) which shows that the fixed point at Uq = U* = is unstable for 2 — d' + 2p > 0. 
Fig. |l| shows a typical RG flow diagram projected onto the (U p ,Uq) plane at d' = 2.5. The 
existence of the off-axis fixed point implies that the nontrivial exponents z and a are given 
by Eq. (11) for d < d' c . 

The analytic formula, Eq. (11) is derived from the one-loop approximation RG trans- 
formation. However, it is in fact exact to all orders. The invariance of the noise spectrum 
D(k) under the scale change, x — > bx, t —>■ b z t, and h —>■ b a h gives the exact exponent 
identity 2p — d! — 2a + z = as argued in 0. This together with the relation a + z = 2 
gives Eq. (11). The difference between the case treated in [[| and in this work is the range 
of applicability of Eq. (11). For the isotropic case, the formula is valid for d! < 3/2 and 
Pmin < P < Pmax with approximately known p max and p m i n . However, for the anisotropic 
case, the formula for the dynamic exponent is valid for d' < d' c = 2 + 2p and p > 0. Hence 
the dynamic exponent z in d' = 2 is nontrivial for finite p. Since the exponents a and z are 
nontrivial in d' = 2, it would be interesting to introduce a stochastic model associated with 
the AKPZ equation having the alternative signs of A's for correlated noise. A possible can- 
didate for such a stochastic model would be a generalization of the 2+1 dimensional Toom 
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model introduced for the white noise || by replacing the random deposition by the spatially 
correlated deposition following the Levy flight distribution. This conjecture is based on the 
analogy to the isotropic case, where the correlated noise was simulated by the correlated 



particle deposition following the Levy flight distribution |K| . The simulations of the model 
is underway O . 



In summary, we have considered the anisotropic KPZ equation with the spatially cor- 
related noise, using the one-loop dynamic RG transformation. In this study, the case that 
the signs of the nonlinear terms are opposite is interesting for both of white and correlated 
noises. Unlike the case of white noise, the stable fixed point exists in the "off-axis" region of 
the parameter space of the effective coupling constants Uq and U p if d! < d' c = 2 + 2p for the 
case of correlated noise. From the stable fixed point, the dynamic exponent z is determined 
as z = 2 — |(p — ^-^) and the roughening exponent a = 2 — z. Therefore, the nontrivial 
values of the exponents a, z are present in d' — 2 for finite p. For dl > d' c , the stable fixed 
point is located at the origin in the parameter space, and the AKPZ equation belongs to 
the Edwards- Wilkinson's universality class. 
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FIGURES 

FIG. 1. Plot of the flow diagram of the RG transformation in the parameter space of U p and 
U for d' = 2.5, p = 0.5 and r A = -2.5. 
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